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The Study of Age-Related Variability of Pigmentation Patterns of the Shells Dreissena polymorpha 
(Bivavia, Dreissenidae) from Different Parts of it's Range. Pavlova, V. V. — Traditionally, shells of 
Dreissena polymorpha Pallas, 1771 mussels are studied as a whole by the investigation of population 
variability of coloration. However, shell surface is divided into zones by the lines of growth delay called 
annual rings, and mussel shell coloration is an array of patterns in the consequent age zones. In the 
present paper traditional term “pattern type” is expanded by its relation to definite shell zone; pattern 
types change along with annual ring formation. The pattern is considered not as a holistic feature but 
rather a sequence of states realized during the mussel’s growth. Frequencies of some pattern types showed 
to be age dependent. Geographic variability of pattern sequences was found. High diversity of pattern 
sequences’ arrays observed in populations may reflect wide adaptation potential of zebra mussels. 
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Visyaenne BO3pacTHOÑ M3MeHYMBOCTM OKpackM pakOBMH MomtocKka Dreissena polymorpha (Bivavia, 
Dreissenidae) u3 pasamaHprx uacteit apeana. [lapnosa B. B. — Tpaguumonno mpm anannse nomyna- 
I{MOHHOTO pasHoo6pasna pucyHKa paKOBMHbI MOJUNOCKOB Dreissena polymorpha Pallas, 1771 paccma- 
TPMBaIOTCA Kak egnHoe genoee. OHAKO NOBePXHOCTb paKOBMH pa3fesAeTCA Ha 30HbI TONOBbIMM KO/Ib- 
yam (IMHMAMN 3af}epxKKM pocta), M OKpacKa CTBOPKM MOJIIOCKa IMpegCTABJIAeT CO6OÑ COBOKYMHOCT 
PMCYHKOB Ha MocmefOBaTebHbIX BOSPAaCTHBIX 30Hax. B FaHHoM paboTe OOujenpMHATOe NOHATME «TMN 
PUcyHKa» JONOMHAeTCA TIpesCTaBeHMeM O IIpMypOueHHOCTH eTO K ONpefeeHHOM BOSpacTHOM 30He, O 
CMeHe TMMOB picyHKa npu oOpasoBaHun rofoBon muHUM. BuepBble pucyHOK paccMaTpMBaeTcA He KAK 
HENOCTHAA XapaKTepUCTUKa, a B BUHE MOCMeOBaTebHOCTM COCTOAHMM, peasIM3yIOLUXCA B Ipolecce po- 
cTa MOJIJIIOCKA. BblaBileHa BO3paCTHad 3ABUCMMOCTb YaCTOT BCTpeyaeMOCTM HEKOTOPbIX THIIOB PUCYHKa. 
Bospacruble MOCMefOBaTebHOCTM TUMOB PUCyHKa NOMBEpPXKEHDI reorpapuyeckow N3MeHYNBOCTN. OT- 
Me¥eHHOE BbICOKOe pasHOobpasie HaOOpoB NOcHeFOBAaTeMbHOCTeM B NOMYIAUMAX MOXKET ABMLATBCA OT- 
PakeHMeM WIMpouaMUx afjalTaljMOHHbIX BO3MOXKHOCTel [[pevicceH. 


Kawuesple Cosa: npeñňccena, OKpacKa paKOBMH, TMI PUCyHKa, BO3pacTHaAA N3MEHYMBOCTÐ. 


Introduction 


Diversity of shell coloration was revealed for many species of both terrestrial and aquatic mollusks. How- 
ever, dreissenids are the only freshwater mussels which have strikingly variable pigmentation patterns. 

Traditionally, description of pigmentation patterns on dreissenid shells includes the discrimination of 
simple elements and their combinations and evaluation of their occurrence frequencies in populations. This 
approach was firstly applied by G. I. Biochino for D. polymorpha (Biochino, 1990). This approach is used by 
A. A. Protasov also (Protasov, Gorpinchuk, 1997). 

G. I. Biochino studied a particular shell as a uniform image. She attributed every shell in a sample to one of 
five coloration types according to prevalent pattern (Biochino, 1990). Opposite to that, A. A. Protasov described 
all characters of shell coloration including background color and elements of pattern which appear individu- 
ally or in various combinations (Protasov, Gorpinchuk, 1997; Protasov, 2000). By this approach the quantity of 
phenotypes in population may reach 80 (Protasov, Gorpinchuk, 1997). 
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Dreissenids have the structures on the outer shell surface, known as “annual rings”, representing the zones 
of growth delay. Visual analysis showed that pattern replacement occurs by the annual ring formation. This 
phenomenon was noted earlier (Marsden et al., 1996; Protasov, Gorpinchuk, 1997), however neither of schemes 
described above took it into consideration (Biochino, 1990; Protasov, Gorpinchuk, 1997). Our approach oper- 
ates the shell not as a whole but as a combination of successive age zones carrying certain pigmentation pattern. 
We suppose that pattern change during the shell’s growth caused by switching of mode of mantle cells that 
are responsible for pigment formation; when studying pigmentation we should take into account these cells’ 
mechanism of functioning. Such idea was applied by S. O. Sergievskiy with co-workers (1995) for investigation 
of pigmentation of gastropod Littorina. Based on this scheme we worked out out’s. 

The aim of the present work is to analyze the pigmentation variability of D. polymorpha from different 
parts of its range taking into account successive pattern change during the mussel growth. 


Material and methods 


Diversity of pigmentation pattern of zebra mussel Dreissena polymorpha from some water bodies in 
Europe and America was studied (table 1). Sampling was carried out by dredge, bottom trawl, rarely by hand. 
All samples were presented by the subspecies D. polymorpha polymorpha except for the sample from Malaya 
Zhemchuzhnaya Banka (northwest part of Caspian), which consisted of specimens of D. p. polymorpha and 
forms intermediate between this subspecies and D. p. andrusovi (Logvinenko, Starobogatov, 1968). Only the 
latter ones were included in the analysis under the designation «D. p. andrusovi». 

Mussels were preserved in ethanol (96 %). In laboratory, they were anatomized, the shells were cleaned 
from tissues, sludge and fouling, and air dried. Only right valves were analyzed since lateral asymmetry of 
pigmentation is typical for zebra mussel (Protasov, Gorpinchuk, 1997). 

Every shell was considered as a sequence of age zones carrying a certain pattern. One zone has one pattern 
type, as there is no pattern change on its area. Thus, in this paper, pattern type is the characters of 
pigmentation of separate age zone of the shell, caused by chromatophore’s work 
in a certain mode that is unchangeable within the whole area of this zone. The part 
of shell of one age carries one pattern type. The older is the mussel the more pattern types its shell has. 

It was shown (Sergievskiy et al., 1995) that mussel pigmentation is determined by disposition and mode 
of work of chromatophore cells and pigments they produce. These investigators described seven modes of 
functioning of chromatophores and the same number of pigmentation variations (table 2). 

Zebra mussel has six types of functioning of chromatophores: 0, 1, 2, 3 a, 3 b, 3 d. Moreover, pigmentation 
diversity is increased as the chromatophore cells are able to operate in different modes in different parts of 
mantle edge, e. g. in the upper region — alternating mode with skew synchronization, in the lower region — 


Table 1. Sampling sites and number of studied mussels 


Ta6anua 1. Mecra c6opa npo6 n konmuecrBo O6paboTaHHOro MarTepnana 





Water body Coordinates Number of mussels 

Volga River, the Rybinsk Reservoir (6 locations)* 58°25' N, 38°20' E 634 
Volga River, the Gorky Reservoir (5 locations)* 57°29' N, 41°21' E 436 
Lake Plescheevo 56°44' N, 38°50' E 51 
The Kama Reservoir 57°56' N, 56°42' E 93 
Lake Forelevoe (near Kaliningrad) 54°43' N, 20°30' E 218 
Northern Dvina River 64°31' N, 40°37' E 77 
The Chograi Reservoir (Kalmykia) 45°34' N, 44°18' E 78 
RR-1 supply channel (Kalmykia) 46°06' N, 41°57' E 48 
Lake Sharony (Kalmykia) 47°53' N, 44°50' E 83 
Ahtuba River 48°41' N, 45°12' E 34 
Northwest Caspian: 

Malaya Zhemchuzhnaya Banka 45°02' N, 48°18'E 240 
Belinskiy Bank 45°48' N, 48°52' E 20 
Danube River (Serbia) 44°50' N, 20°31' E 124 
The Perućica Reservoir (Serbia) 43°58' N, 19°21' E 186 
Lake Frie (USA) 41°48' N, 82°59' W 43 
Lake Michigan (USA) 43°11' N, 86°26' W 98 
Total 2463 








* Coordinates of center of water body are given. 
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Table 2. Mode of chromatophore cells work according to parametric system of mollusk shell pigmenta- 
tion description (Sergievskiy et al., 1995) and designations of pattern types in zebra mussel (after Biochino, 
1990 and Protasov, Gorpinchuk, 1997) 


Ta6anua 2. Pexumbl þpyukynonnposanna xpoMaToopos cormacHo napamerpnueckoÑñ CucTeMe onn- 
CaHMA OKpackM paKkOBMH MOMIIOCKOB (CepriueBckuii n gp., 1995) n cooTBeTCTByIoMIMe MM O6osHaTeHHA 
THMOB (WIN I1eMeHTOB) PUCyHKa pakOBMH gpeřccen (mo cxemam: Buounno, 1990 n IIporacos, Topmmy- 
uyk, 1997) 














Mode of chromatophore cells functioning Designations | Designations (after 
(after Sergievskiy et al., 1995) Shell pigmentation | (after Bio- Protasov, Gorpin- 
chino, 1990) chuk, 1997) 
1) Functioning of all chromatophores in Dark coloration OO D, 
one constant mode — monochrome shell without pattern 
coloration 
2) Functioning of some of chromatophores Longitudinal stripes RR LE 
in one constant mode 
3) Functioningof a) unstable Irregular spots — L 
chromatophores in synchronization; 
alternating mode p) single-phase Transversal stripes AA G 


with periodic 


; synchronization; 
synthesis delay. . “ P : : 
Different types of c) antiphase Chessboard Not registered in zebra mussel 
synchronization of synchronization; 
pigment cell work: d) skew Broken/zigzag lines CC F,H,J,K 
synchronization. 
Absence of chromatophores or their Depigmentation DD C, 


inactivity (this mode was not designated 
by Sergievskiy with co-authors; we use 
character «0»). 





alternating mode with single-phase synchronization. Traditionally, such complex pattern of dreissenid mussels 
is described dividing it into simple elements. We used the same approach. 

Pigmentation pattern of dreissenid shells is composed of transversal and longitudinal elements that are 
designated as follows: 

Transversal elements (we use designations proposed by G. I. Biochino (1990), with additions): 

A — regular arcuate pattern along the whole shell height (3 b type by Sergievskiy et al.,1995); 

C — zigzag pattern along the whole shell height (3 d); 

AC — intermediate pattern between A and C: the upper shell half carries zigzag, the lower — arcs 

(combination of 3 b and 3 d); variations of this element: 

ACs — regular arcs along the nearly whole shell height, curves are only on the uppermost part of the 
shell; 

ACi — regular arcs along the nearly whole shell height, curves are only on the lowest part of the shell; 
Ad — arcs with teeth in the middle; 

D — light shell without pattern; this type was put to transversal elements conditionally (0 type). 

Longitudinal elements: 

R — radial lines (2 type); variations of this element: 

Rd — single light radial line; 

Ro — single dark radial line; 

Rm — multiple thin radial lines, light and dark alternately; 

F — border between different elements running longitudinally. 

Pattern type is an ensemble of elements and their enumeration is the designation of pattern type. 

For the description of pattern type sequences, age of mussels was determined by annual rings and mussel 
length distribution in a sample. Sequence of pattern types from age 0+ was revealed (fig. 1). After 3+ age change 
of pattern types was not registered, therefore the record was made till this age. 

Data from individual shells were summarized. Frequencies of pattern types in age zones 0+, 1+, 2+, 3+ 
were obtained for every sample. Term “age zone ...+” designates shell parts of mussels of different ages formed 
at this age, e. g., phrase “age zone 0+” means shells of underyearlings (0+) and shell parts before the first annual 
ring (age of 1 year) of older mussels. “Age zone 1+” means shell area after the first annual ring in yearlings (1+) 
and shell parts between the first and the second annual rings (age of 2 years) of older mussels etc. 

Frequencies of different pattern types were determined both in different age zones separately and in to- 
tal, summarizing data for all zones. Spearman rank correlation coefficient was used for the analysis of rela- 
tion between number of pattern types and geographic latitude and between pattern type frequencies and age. 
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3+ 2+ I+ 0+ 
A A AC C 


Fig. 1. Change of pattern types on zebra mussel shell. 
The present shell has four age zones (0+, 1+, 2+, 3+). The pattern sequence is C-AC-A-A. 


Puc. 1. Cmena THN0B pucyHKa Ha pakosnne D. polymorpha. 
Jannaa CTBOpKa MMeeT YeTbIpe BO3PacTHbIe 30HbI (0+, 1+, 2+, 3+). 3ammcb MOcMefOBaTebHOCTM TUMOB MMeeT 
Bng C-AC-A-A. 


Similarity between populations was estimated following Zhivotovskiy (1982) (r ,,,, index). The significance of 
differences between populations was determined using y’ criterion (Lakin, 1980). Correlation between popula- 
tion similarity and geographic distances between them were estimated using Mantel test: GenAIEx 6 (Peakall, 
Smouse, 2006); significance level was determined using 9999 permutations. 

Diversity of sequences of pattern types was estimated for whole samples, for mussels of all ages. For every 
population share of individuals carrying different sequences was calculated. 


Results 


Change of pattern type during mussel growth is typical for zebra mussel but no pattern 
change was revealed in 3.9 % of 2+ mussels and 0.9 % of 3+ mussels. Pattern types in adja- 
cent age zones may be either the same or different (e.g. in shell with sequence AC-AC-A-A 
pattern remains the same by switch from 0+ to 1+ and from 2+ to 3+ and changes by switch 
from 1+ to 2+). The share of specimens with either constant or changing patterns varies 
between samples. On average pattern remains the same in 16.6 % of individuals by switch 
from 0+ to 1+, in 32.5 % of individuals by switch from 1+ to 2+, in 65.9 % of individuals 
by switch from 2+ to 3+. Thus, pattern types in older age zones tend to be more stable than 
in the younger. 


Total frequencies of pattern types 


In total, 16 pattern types were found in zebra mussel (fig. 2, table 3). 

Pattern types C, AC, Ad, A are the main types present in all populations (all other types 
are secondary). Total frequency of main types varied from 73.8 to 100 % across different 
samples. Pattern was presented only by basic types in 6 from 26 samples (4 locations in 
the Rybinsk Reservoir, in Ahtuba River and in Belinskiy Bank). The maximum diversity 
(8-9 pattern types) in D. p. polymorpha was found in samples from the Kama and Chograi 
Reservoirs, from the Danube River and Lake Michigan. Very high variability was obtained 
in «D. p. andrusovi» from Malaya Zhemchuzhnaya Banka (MZhB) (14 pattern types). 

Average within-sample pattern diversity in D. p. polymorpha is 6 pattern types. There 
is no significant correlation between number of pattern types in the sample and geographic 
latitude (R = -0.37, p > 0.05) (fig. 3). 

Mean intersample similarity r ,,.. is 0.892 + 0.05. Minimal value was found in pair 
Danube River — Lake Michigan (0.441), maximal value — for two samples from the Gor- 
kiy Reservoir (0.996). Correlation between r „ and geographic distance between locations 
is statistically significant (R,, = -0.45; p = 0.002). 
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Fig. 2. Pattern types on zebra mussel shells. 


Puc. 2. Tunai pucyHka Ha pakosnnax D. polymorpha. 
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Table 3. Frequency of occurrence of pattern types in zebra mussel populations 
Ta6anuya 3. Bcrpeaaemocrb TUNOB pucyHka B Honynanuax D. polymorpha 
Pattern type Frequency, N (%) of samples Share in different samples, min-max, % 
C 26 (100.0) 5.5-58.5 
AC 26 (100.0) 5.2-50.0 
Ad 26 (100.0) 0.5-37.6 
A 26 (100.0) 0.5-67.1 
ACs 16 (61.5) 0.2-5.6 
ACi 13 (50.0) 0.4-5.2 
D 8 (30.8) 0.5-9.1 
AFD 1 (3.9) 13 
ARd 4 (15.4) 0.5-2.5 
CRd 5 (19.2) 0.3-2.4 
ACRd 2 (7.7) 0.5-0.9 
ARm 4 (15.4) 0.6-9.4 
ACRo 1 (3.9) 0.2 
DRo 1 (3.9) 21 
ARo 3 (11.5) 1.7-1.9 
ARoRm 3 (11.5) 0.2-3.5 
10 
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Fig. 3. Latitudinal variability of the number of pattern types. Samples: 1-6 — the Rybinsk Reservoir; 7-10 — 
the Gorky Reservoir; 11 — Lake Plescheevo; 12 — the Kama Reservoir; 13 — Lake Forelevoe; 14 — Northern 
Dvina River; 15 — the Chograi Reservoir; 16 — RR-1 channel; 17 — Lake Sharony; 18 — Ahtuba River; 19 — 
Belinskiy Bank; 20 — Danube River; 21, 22 — the Perućica Reservoir; 23 — Lake Erie; 24, 25 — Lake Michigan. 


Puc. 3. Msmenenne komMuecTBa THNOB puicyHKa c reorpapuyeckon mnporoň. O6osHayenna BEIOopoK: 1-6 — 
Ppi6uucKoe ByXp.; 7-10 — ToppKoscKoe Byxp.; 11 — o3. Inenyeego; 12 — Kamcxoe Baxp.; 13 — 03. Popene- 
Boe; 14 — p. Cepepuaa {Bnna; 15 — Horpañckoe Byxp.; 16 — Kanan PP-1; 17 — os. Iapouni; 18 — p. Axty6a; 
19 — Benuucxun 6anx; 20 — p. HyHan; 21, 22 — syxp. Ilepyauma; 23 — 03. Əpn; 24, 25 — 03. Mvanran. 
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Fig. 4. Frequences of main pattern types at different age zones on zebra mussel shells. 


Puc. 4. Uacrorbi OCHOBHDBIX TMMOB piicyHKa Ha pa3/IM4HbIX BO3PaCTHbIX 30HaX pakoBMH D. polymorpha. 


Frequencies of pattern types on different age zones 


General regularities 

The highest pattern diversity was revealed in early age zones, the lowest — in 3+ zone. 
In 0+ age zone 11 pattern types were found in D. p. polymorpha and 12 in «D. p. andrusovi»; 
in 1+ zone 14 and 10 pattern types, respectively; in 2+ and 3+ zones — 11 and 6 pattern 
types in D. p. polymorpha («D. p. andrusovi» was presented only by young mussels). 

Some pattern types have been shown to be age-related. Type C is specific for early age 
zones, only one individual had this type in 3+ zone (fig. 4). Correlation between frequency 
of this type and the age is 0.83 (p < 0.05). The frequency of type A rises as age increases (R = 
= 0.83, p < 0.05). Age-dependence of frequencies of types AC and Ad is weaker (R is 0.25 
and 0.52, respectively, p < 0.05). 

Regularity of age alteration of main types C, AC, Ad, A described above is basic in 
D. polymorpha and is obtained in all samples (fig. 5). 

Mean intersample similarity of pattern on separate age zones decreases from 0+ to 3+ 
(T piv 18 0.856 + 0.006, 0.856 + 0.006, 0.780 + 0.012, 0.730 + 0.020, respectively). 


Local characteristics 

Figure 5, a shows the most typical pattern types distribution on age zones. It is typical 
for all samples of D. p. polymorpha from the Rybinsk and the Gorkiy Reservoirs and from 
Lake Forelevoe, the Perućica Reservoir and the Danube River. Frequencies of secondary 
types did not exceed 8.9 % in these samples. 

Population from Lake Plescheevo is evidently different from the Upper Volga basin 
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samples despite of it is the source of penetration of zebra mussel in the lake (Bakanov, 
1992). The most prominent difference is the frequency of type Ad on age zones from 0+ 
to 2+ which is the highest among all samples (it achieved 51.1 %). Occurrence of type C 
on age zone 0+ is one of the lowest. Similarity of populations from Lake Plescheevo on the 
one hand and from the Rybinsk and the Gorkiy Reservoirs on the other in 0+ varies from 
0.799 no 0.934 (differences between sample from Lake Plescheevo and any sample from 
two reservoirs are statistically significant; x? is 7.5-37.7, number of degrees of freedom is 
1-2, p < 0.05). However, the similarity of mussel samples from different locations in two 
reservoirs is 0.894-1.000. Upper Volga populations show more likeness to geographically 
distant populations from Lake Forelevoe and the Kama Reservoir (r „is 0.939-0.994 and 
0.881-0.950, correspondingly) than to populations from nearby Lake Plescheevo. 

Sample from the Northern Dvina River is the only where type C is prevalent not only 
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Fig. 5. Examples of fluctuation of pattern type frequencies in zebra mussel populations. 


Puc. 5. Pacnpegenenne uacToT TUNOB puncynKa B Honynaynax D. polymorpha. 
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in 0+, but also in 1+ (85.7 and 76.3 %, correspondingly) (fig. 5, b). 

High frequency of type A in early age zones is typical for many populations (Belinskiy 
Bank, MZhB, RR-1 channel, the Chograi and Kama Reservoirs, Lakes Erie and Michigan). 
It rises 40 % in 0+ (Lake Erie), in 1+ varies from 31.1 to 63.0 %, on older shell zones is 
72.1-100 %. In the samples mentioned low frequency of type AC is characteristic as well; it 
does not exceed 33.3 % at different ages. The portion of secondary types is high, it reaches 
18-20 % (in samples from RR-1 channel and the Chograi Reservoir). In the population of 
«D. p. andrusovi» (MZhB) total frequency of secondary types is 31.9 % in 0+ and 20.2 % in 
1+. Unique types ACRo and DRo were found in this population. 

Specific feature of samples from Lake Sharony and the Ahtuba River is a high share of 
type AC. Unlike other samples, this type is frequent in not only 0+ and 1+ age zones (up to 
56.1 and 60.0 %), but also in 2+ and 3+ (up to 66.7 u 50.0 %). 

Thus, the differences between zebra mussel populations are manifested in changes of 
frequencies of main types C, AC, Ad and A with the age. The secondary types show weak 
age dependence. 


Age-related sequences of pattern change 


Pattern types substitute each other during the zebra mussel growth. Change of pattern 
types is not chaotic. Usually the pattern changes from zigzag to arched as mussel grows. 
Figuratively, lines of pattern “are straightened” (fig. 1). As the sequence of types it is encod- 
ed as C-AC-Ad-A. Some stages in this row may be absent, while others may be repeated, 
but the general order of types change from zigzag to arch is always present. Inverted order 
is extremely rare and may be considered an aberration. The pattern type A may be consid- 
ered as “final”: this pattern does not change after its realization. Some shells carrying this 
type from 0+ and during the whole ontogenesis were found. 

The diversity of pattern types determines the variety of sequences in population. For 
example, the greatest variety of sequences (25 at age 1+) was revealed in the sample of 
“D. p. andrusovi’, in which the number of pattern types is the greatest. In D. p. polymorpha 
samples number of sequences do not exceed 13; it is 6-8 in those ones where only main 
pattern types were revealed. 

The array of sequences in population consists of a big amount of rows among which 
one or two are dominant (fig. 6). 

Comparison of sequences arrays has shown that the most similar were populations 
from the Rybinsk, Gorkiy, Kama and Perućica Reservoirs, from Lake Forelevoe and Danube 
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Fig. 6. Scheme of pattern change in zebra mussel population (by the example of the Perućica Reservoir). Data for mussels 
of ages from 1+ to 3+ are aggregated. 

Numbers denote the percentage of individuals possessing the given sequence at the given age. Dotted lines 
denote the pattern change not revealed in the mussels of current age but revealed in older ones. 


Puc. 6. IlocneqopatembHOcTM CMeHbI pucyHKa B Monymayun D. polymorpha (ua mpumepe Byxp. Ilepyumua). 
O6pequHéHHaa CxeMa [JIA MOJIIOCKOB Bo3spacToB l+ - 3+. 

Lndpamu 0603Ha4eHbI YacTOTbI BCTpedaemocTu Ocobert (B IpoleHTax), y KOTOPbIX OTMeyeHa JaHHaa 
TIOCMe{OBATEIbHOCTb B JAHHOM BO3pacTe. IlyHKTMpoM OOO03HAaYeHDI HepexObl, He OTMeYCHHbIe B TeKYLIEM 
BO3pacTe, HO BbIABJICHHbIe Y MOJIIKOCKOB CTApHIMX BO3pacToB. 
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Table 4. Sequences of pattern types for zebra mussel populations 
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Ta6anua 4. IlocneyoparenbHocrn CMeHbI THNOB pucynkKa B nonyrarnax D. polymorpha 














Water body Main sequences Secondary sequences 
The Rybinsk Reservoir C-AC-Ad C-C-AC-A 
C-AC-A-A AC-AC-Ad-A 
—~A-A- 
The Gorky Reservoir C-AC-Ad AC-A-A 
C-AC-A-A C-Ad-A 
AC-Ad-A 
The Kama Reservoir -«- C-Ad-A 
Lake Forelevoe -«- AC-Ad-Ad 
AC-A-A 
C-AC-AC-Ad 
Northern Dvina River C-C-A-A C-C-Ad-A 
C-C-Ad-Ad 
Lake Sharony C-AC AC-Ad 
AC-A C-A 
AC-AC 
The Chograi Reservoir C-Ad AC-A 
C-AC Ad-A 
C-A 
RR-1 channel C-A AC-Ad 
C-AC-A 
Ahtuba River C-A AC-AC-AC 
C-AC-AC 
Malaya Zhemchuzhnaya Banka, Ad-A AC-A 
“D. p. andrusovi” D-ARm C-A 
C-A ARd-A 
AC-A DRo-A 
ARo-AC 
Belinskiy Bank C-A C-AC-A 
C-Ad-A 
Danube River, the Perucica Reser- C-AC C-Ad 
voir C-AC-Ad-A 
C-AC-A-A 
Lake Erie ?-A-A-A 
Lake Michigan ?-Ad-A-A 





Note. «?» means damaged shell zones where no pattern can be recognized. 


River. They are different only in secondary frequencies (table 4). 

The specific feature of the Northern Dvina population is the predominance of C-C 
sequence when 0+ is changed by 1+. 

Pattern types change C-A is typical by the switch from 0+ to 1+ for south water bodies. 
The highest frequency of this sequence is found in samples from the Ahtuba River and Belin- 
skiy Bank. It was also noted in MZhB, Chograi Reservoir, Lake Sharony and RR-1 channel. 

The set of sequences in «D. p. andrusovi» is specific, it includes unique rows D-Arm, 
ARd-A, DRo-A and ARo-AC. 

Main sequences of pattern types in North American samples of zebra mussel were 
failed to be obtained since upper pigmented layer of the shells was eroded on anterior shell 
part in many individuals. Nevertheless it was noted that the most frequent sequences be- 
ginning from 1+ age zone are A-A-A and Ad-A-A, and that is similar to populations from 
the South Volga basin. 

Therefore, fast “straightening” of pattern lines (from C to A in 1+ age zone, without 
intermediate stages) occurs in southern zebra mussel populations from the Volga River 
basin; “straightening” is slower in northern ones (through intermediate types AC and Ad 
in 1+ - 2+ age zones); in the most northern sample (the Northern Dvina River) type C is 
obtained during the longest period and “straightening” begins from 2+ zone. 
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Discussion 


Mussel shell pigmentation is determined by the operation of mantle chromatophore 
cells. Pattern appears as a result of functioning of these cells in a certain mode (Sergievskiy 
et al., 1995; Boettiger et al., 2009). Factors that induce growth delay (slowdown of shell 
secretion by mantle cells) may also affect pigmentation-responsible cells altering mode of 
their functioning. Such factors are reduced or elevated temperature, spawning and others 
(L’vova, 1980; Allen et al., 1999; Karatayev et al., 2006). 

Obviously, chromatophores’ functioning is influenced by environmental conditions. 
Interpopulation differences in frequencies of pattern types and their sequences found in 
D. polymorpha (fig. 5, table 3) confirm it. Apparently, degree of population similarity is also 
determined by the intensity of gene flow between them. For example, differences in pattern 
types frequencies between Upper Volga reservoirs are low since constant gene exchange 
occurring between them through veligers. Population from nearby Lake Plescheevo 
possesses dissimilar pattern types frequencies, that is caused by isolation (there is no stable 
interconnection between Upper Volga populations and one from Lake Plescheevo). 

It was noted earlier that pattern variability in zebra mussel is geographic-related (Biochino, 
1994; Protasov, 2000). It was shown that frequency of arched pattern (type A) in southern 
population in Europe is higher than in northern. Our results support this and specify that it 
is happening due to appearance of type A on younger shell zones in southern populations. 
Consequently, type A is found in a larger number of mussels as it is final type. It is assumed 
that phenotypes are differentially adaptive to the environmental conditions (Yablokov, Larina, 
1985). For example, Protasov with co-authors (1997) revealed that individuals of D. polymorpha 
with arched pattern were more resistant to hypoxia, than specimens with zigzag pattern. 

Described age-relation of pattern variability in dreissenid shells allows assuming that the 
highest adaptability to several factors develops successively during the ontogenesis. Observed 
diversity reflects enormous adaptation resources of dreissenid populations: an individual is 
best adapted to a certain factor at a certain age. As a result, the population possesses adaptive 
potentials to wide range of conditions, even if they change constantly. Pattern diversity may serve 
as explanation for high adaptability of zebra mussel and its success as invasive species. 

Despite of variability of the pattern types sequences, “straightening” of pattern lines is 
specific for most of mussels. The most curved lines (type C) are typical for mussels at early 
stages; older age zones of most of individuals carry regular arched pattern lines (type A). 
Probably, skew synchronization of chromatophore cells work (Sergievskiy et al., 1995) 
forming type C is connected with intense mussel growth. Growth rate slows with age, and 
possibly stabilization of chromatophores’ work in a single-phase mode occurs when all cells 
of mantle edge switch simultaneously and arched pattern occurs. 

This general regularity does not contradict the conclusion mentioned above on the 
variety of adaptive possibilities of populations since deviations exist. First of all, initial 
pattern type is not always C. Pattern type variability in age 0+ is one of the highest in 
populations. In some cases the initial type is A and it is also “final”, and hence persists 
during the whole ontogenesis. Moreover, succession of transition from C to A differs in 
populations (table 3). 


Conclusion 


Totally, 16 pattern types were found in zebra mussel. The pattern variety in intermediate 
forms between D. p. polymorpha and D. p. andrusovi (“D. p. andrusovi”) is more than that 
in D. p. polymorpha. General regularity of age-related pattern change was revealed: decrease 
of type C portion and increase of type A portion during the transition from early to older 
age zones. The rate of this “straightening” grows from north to south. 


Author appreciates Yu. V. Gerasimov, I. A. Stolbunov, A. N. Kasyanov, D. D. Pavlov, D. E. Vekhov, 
E. V. Nikitenko, D. A. Guseva, I. V. Pozdeev, T. F. Nalepa, B. Mickovic for help in mussel sampling. 
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